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A few years ago, Campos investigated the critical phenomena of thick branes in warped spacetimes 
[Phys. Rev. Lett. 88 (2002) 141602]. Inspired by his work, we consider a toy model of thick branes 
generated by a real scalar field with the potential V{4>) = a(j)^ — h<j)^ + c0®, and investigate the 
variation of the mass parameter a on the branes as well as the localization and resonances of 
fermions. An interesting result is found; there is a critical value for the mass parameter a, and 
when the critical value of a is reached the solution of the background scalar field is not unique and 
has the shape of a double kink. This happens in both cases with and without gravity. It is also 
shown that the numbers of the bound Kaluza-Klein modes of fermions on the gravity-free brane and 
the resonant states of fermions on the brane with gravity increase with the value of a. 



I. INTRODUCTION 

The idea of embedding our Universe in a higher di- 
mensional space-time provides new insights for solving 
some long-standing problems in physics, such as the prob- 
lems of the gauge hierarchy [l--!] and cosmological con- 
stant [4]. Since the discovery of four-dimensional grav- 
ity can be realized on a brane in fivc-dimcnsional (5D) 
space-time [•'»], the brane- viforld theories have received 
considerable attention. The branes in Randall and Sun- 
drum model [3, 5] are fixed on some points along the 
fifth dimension, -with a form of 5 function. The Ran- 
dall and Sundrum brane--world model is very ideal, it 
is infinitely thin, and the formation of it lacks the dy- 
namical mechanism. A natural extension of thin brane 
models is the thick brane scenarios with the gravity cou- 
pled to five-dimensional scalar fields . Furthermore, 
there are some other thick brane models such as the mod- 
els -with thick branes arising from pure geometry [34-39] 
or fermion self-interaction [4IJ, 41]. Some reviews about 
thick brane- viforld models can be found in Refs. [42-45]. 
In brane--world theories, possessing the possibility of lo- 
calizing all kinds of matter fields on branes is very im- 
portant. The localization of particles on a thick brane 
-without gravity -was suggested a long time ago [4] . In the 
case -with gravity, the localization of fermions and bosons 
on various brane -worlds has been studied in much of the 
literature [4()-G9]. Furthermore, for some thick branes, 
the massive fermions trapped on them are unstable and 
have a finite confinement lifetime [53, 59, 68, 70-74]. 

Recently, Campos discussed the critical phenomena of 
thick branes in -warped spacetimes [14]. In his -work, the 
thick branes are constructed with a complex scalar field, 
and the real part of the complex scalar field is used as 
an order parameter field to describe the first-order phase 
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transition. In his paper the potential of the background 
complex scalar field $ has the form of 



(1) 



where 0^ and (pi are the real and the imaginary parts of 
the complex scalar field $, respectively. Taking the effect 
of the temperature into account, the mass parameter A 
should be a function of the temperature T [75-77]. Cam- 
pos found that the presence of gravity would lower the 
critical temperature of the phase transition correspond- 
ing to what occurred in Euclidean space [14]. And when 
the critical temperature is reached an interface that inter- 
polates two bulk ordered phases splits into two interfaces 
with a disordered phase between them [14]. 

Inspired by the work of Campos, in this paper, we con- 
sider the thick branes with a real scalar field as the order 
parameter and the background field. For the purpose 
of investigating the universal property of the model, the 
potential of the scalar field is taken as 



(2) 



where a,b,c > 0. Taking the effect of the temperature 
into account, different from the 0* model, then a and b 
should all be functions of temperature T [78, 79]. But to 
calculate the exact forms of a{T) and b{T) is very diffi- 
cult in the five-dimensional curved space-time, although 
a fivc-dimcnsional field in the bulk can be regarded as 
a four-dimensional field on the brane with an infinite 
Kaluza-Klein (KK) tower of masses [80, 81]. In brane- 
world theories, we have not found the exact form of the 
effective potential of the (j>^ model with temperature cor- 
rections. 

Without the exact forms of a{T) and b{T), here we 
just investigate the effects of the variation of the mass 
parameter a on the fermion localization and resonances. 
This can partly represent the effects of the temperature 
on the branes in this (p^ model. In our investigation, 
we find that the number of the bound states or of the 
resonant states increases with the parameter a. Besides, 
when the critical value of a is reached the solution of 
the background scalar field is not unique and it has the 
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shape of a double kink, which happens in the two cases 
with and without gravity. That means when the system 
approaches the critical value of a the width of the brane 
will become uncertain. 

The paper is organized as follows. In Sec. II, the effects 
of the variation of the mass parameter a on thick brancs 
are discussed in the case without gravity. This includes 
the effects on the solution of the background scalar field 
and on fermion localization. In order to localize fermions 
on the branes, the Yukawa coupling between the back- 
ground scalar field and fermions is introduced. Then 
Sec. Ill deals with the case of thick branes with grav- 
ity. For the same Yukawa coupling, the gravity on the 
branes will lead to the disappearance of the discrete spec- 
trum of fermions, and only one bound state of fermions, 
i.e., the four-dimensional massless fermion, resides on the 
branes. Furthermore, fermion resonances with a finite 
lifetime could appear on the branes, and their number 
increases with the mass parameter a. Finally, a brief 
conclusion is presented in Sec. IV. 



II. WITHOUT GRAVITY 

We begin with the case without gravity. Consider a 
5D model containing one real scalar field. The action is 



d xdy 



1 



(3) 



where M,iV = (0,1,2,3,5) and 77" 
diag(— 1, 1, 1, 1, 1). The self-interaction potential of 
the background scalar field is of the form as (2). There 
are three local minima for V{(j)), one is at 0'^^ = 
corresponding to a disordered bulk phase and the other 
two are at (/)(^^ = — (/)('^) = v with 



\/\P- — 3ac 
3c 



3c ' 



(4) 



they are degenerated and correspond to ordered bulk 
phases. When a = Oc [a^ = 6V4c], V{<^^^^) = 1^(0^^') = 
has three degenerated global minima. So 
we take flc as the parameter corresponding to the critical 
value of a [14]. The profiles of the potential as a func- 
tion of the are shown in Fig. 1. This is very similar to 
the model in Ref. [14], where the thick branes are con- 
structed from a complex scalar field and the action has 
Z3 symmetry. The usual hypothesis is that is only the 
function of the extra dimension y, i.e., 4> = 4>^V)- 



A. Effects of the variation of the mass parameter a 
on the thick brane solutions 



The equation of motion for </> corresponding to the ac- 
tion (3) reads as 



dV{^) 



(5) 




FIG. 1. The profiles of the potential as a function of 4> 
with h = 2 and c = 1 (oc = 1). The lines correspond to 
a — 0.8, 0.837, 0.9, 1 from bottom to top, respectively. 



where the prime stands for the derivative with respect 
to y. When y — >• ±00, (j){±oo) = ±v and 0'(±oo) = 0. 
Considering the symmetry of the potential F(</>), we can 
set 0(0) = 0. In summary, the boundary conditions of 
Eq. (5) are 



(+00) = -0(-oo) = V and 0(0) = 0. 



(6) 



We have numerically solved the equation (5) using the 
above boundary conditions with 6 = 2, c = 1 and differ- 
ent values of a with a < = 1; the results are shown in 
Fig. 2. From the figure we can find that with the increase 
of the value of a the profile of the background field has 
the tendency to hold the form of a double kink. For the 
energy density of the background scalar field, we have 



1 



P = o0 + 1^(0)- 



(7) 



Its profiles with different values of a are shown in Fig. 3. 
We sec that for a = 0, the energy density has a single- 
peak around y — 0. As the value of a is increasing, 
the maximum of p splits into two new maxima, and the 
distance of the new maxima increases with the value of 
a. That means the brane gets thicker and thicker with 
the increase of the value of a. 

When a = Oc we find that the profile of has a 
double-kink form and the solution of is not unique [8 2], 
which means that the boundary conditions (6) can not 
determine the solution uniquely. Figure 4 shows the 
profiles of the different solutions of with a ~ ac- 
This is very different from the noncritical case. So in 
the brane-world model, when the critical value of a is 
reached the ordered-ordered brane (located between two 
ordered phases) splits into two ordered-disordered inter- 
faces (with each one located between an ordered phase 
and a disordered phase). The phase transition can be 
identified by 0'(O) = 0. 
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FIG. 2. The profiles of tfie bacfcground scalar field (f) in the 
gravity- free case with 6 = 2 and c = 1 (ttc = 1). The three 
lines correspond to a = 0, a = 0.9, and a = 0.99, respectively. 
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FIG. 3. The profiles of the energy density for the background 
scalar field in the gravity-free case with 6 = 2 and c = 1 
(flc = 1). The three lines correspond to a = 0, a = 0.9, and 
a = 0.99, respectively. 




FIG. 4. The profiles of the background scalar field (f) in the 
gravity- free case with a = Uc- The three lines correspond to 
three numerical solutions with different initial guess solutions. 



B. Fermions on the brane without gravity 

In order to localize fermions on the brane, we need to 
introduce the Yukawa coupling of the fermion with the 
background scalar field (j). The action of the bulk fermion 
is 



(8) 



where rj is the coupling constant, 7^^ = (7'', 7^) {n = 
0, 1, 2, 3) with 7'^ the usual 4D Dirac matrices and 7^ the 
chirality operator. From the above action we can obtain 
the equation of motion for the fermion 



* = 0. 



(9) 



In order to solve this equation, we need to separate vari- 
ables with the usual KK and chiral decomposition: 



E 



{i'l.n[x),fl,n{v) +i'Rn{x)fRn{y)y (10) 



where the sum also includes the integral over all contin- 
uum cases. The 4D left- and right-handed fermions iphn 
and ipRn satisfy the following coupled equations: 



The KK modes and /r,i satisfy 

fhn + V^f^f Ln = ninfRn, 

and the following orthonormality conditions: 



(11) 
(12) 

(13) 
(14) 



fhmlRndz 0. 



fRmfRndz, 



(15) 



From Eqs. (13) and (14) we can get the following 
Schrodinger equations with eigenvalue m^: 

[-a2 + Vi(y)]/L„ = m2/L„, 

[^d^y+VKiy)]fRr. 

where the potentials are 



mlfRn, 



(16a) 
(16b) 

(17) 
(18) 



Before solving the Schrodinger Eqs. (16), let us take a 
look at the behavior of Vl(2/) and Vji{y) at y = and 
y — >■ ±00. 

We first focus on the case of a < Oc. With the bound- 
ary condition ^(0) = 0, we have Vl(0) = —rj<j)'{Q) and 
^r(O) = +ri4>'{0). From the numeric solutions in Fig. 2, 
we can find that 0'(O) > 0. So in the case of 77 > 0, 
Vl(0) < and Vr,(0) > 0. Without loss of generality 
we choose 77 > 0. When y — ±00, c/)(±oo) = ±1; and 
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a 


The maximum of Vl 


'2 

mo 


'2 

nil 


2 

ma 


'2 





5.333333 





4.797808 






0.5 


4.774852 





3.748817 






0.9 


4.1869012 





2.118862 


4.003805 




0.99 


4.019852 





0.995953 


2.496547 


3.824792 



TABLE I. The mass spectrum of the left-handed fermion in 
the gravity- free case with different values of a and fixed b — 2, 
c — 1, and rj — 2. 



a 


The maximum of Vr 




mi 







5.333333 


4.797808 






0.5 


4.774852 


3.748817 






0.9 


4.186901 


2.118856 


4.003659 




0.99 


4.019852 


0.994289 


2.496545 


3.824551 



TABLE II. The mass spectrum of the right-handed fermion 
in the gravity-free case with different values of a and fixed 
& = 2, c = 1, and r) = 2. 



III. WITH GRAVITY 

In the scenario of the thick brane world with gravity, 
the action is as fohows: 

S = j d^xdy^g ( - A - \y'dM<l>dN4> - ^(0)) , 

(19) 

where A is the 5D bulk cosmological constant, and the 
scalar potential V{4>) is the same as (2). The metric in 
this model is assumed as 

ds^ = e^^'^y^ri^^dx^'dx'' + dy^, (20) 

where fj,,^ ~ 0, 1, 2, 3, rj^j^i, ~ diag(— 1, 1, 1, 1), and e^"*^^^ 
is the warp factor. The usual hypothesis is that A is only 
the function of the extra dimension y. 

A. Effects of the variation of the mass parameter 
a on the thick brane solutions 



(/)'(±cxd) = 0, so Vl,r(±oo) = ?7^w^, which is a constant 
greater than zero. 

The shapes of the potentials Vl and Vr for different 
values of a are shown in Fig. 5. From the figure, we can 
find: (a) The potential of the left-handed fermion has 
a negative potential well, so the left fermion may has a 
massless mode, (b) As a is close to the critical value 
there will be two minima in the potential well, (c) For 
the right-handed fermion case, the potential Vr^ has also 
a finite well when rj > (j)'(0)/v and does not have any 
potential well when < 77 < (j)'{0)/v, but the minimum 
of the potential is always greater than zero for any 77 > 0, 
so there will not exist a zero mode. 

Here, it should be noted that the form of the potential 
well is also related to other parameters such as b, c, and 
77. But in this paper we lay the emphasis on the effects of 
the variation of mass parameter a on the branes, so the 
effects of other parameters will not be taken into account. 

We have solved Eqs. (16) numerically for the case 
of the value of a below its critical value. And we get 
the mass spectra of the left- and right-handed fermions, 
which are shown in Tables I and 11, respectively. From 
the two tables we find: (a) The number of the bound 
KK modes of both the left- and right-handed fermions in- 
creases with the value of a. (b) The spectra of the massive 
KK modes of the left- and right-handed fermions are al- 
most the same, which demonstrates that a Dirac fermion 
could be composed from the left- and right-handed bound 
KK modes [72]. 

For the case of a = Oc, because the solution of the 
background scalar field </> is not unique and there are no 
other conditions to fix it, so the localization of fermions 
in this case is not taken into account here. 



The equations generated from the action (19) with 
ansatz (20) reduce to the following coupled nonlinear dif- 
ferential equations 

0" = -W-^^, (21) 

A" = -^r, (22) 

A'^^h^'^-2V)-^. (23) 
6 3 

Equation (23) can be used to fix the bulk cosmological 
constant A. When y — > ±00, |(/)(±(X))| — cj)'--^^ = — 0^'^' = 
V. In this paper we chose <^(-|-oo) = v and (t>(—oo) — —v. 
The boundary conditions can be read as follows [14]: 

A{0) = ^'(0) = 0(0) = 0, (24a) 
0(+oo) = -0(-oo) = V. (24b) 

With the boundary conditions (24) we have solved Eqs. 
(21) and (22) numerically. The results arc shown in 
Figs. 6 and 7. And the profiles of the matter energy 
density 

Too^e'^'ilr' + Vm (25) 

are shown in Fig. 8. The effect of the variation of the 
mass parameter a on the energy density is similar to the 
case without gravity. Through the comparison of Figs. 3 
and 8 we see that the existence of gravity can change the 
behavior of Too a-t infinity. When y — > ±00, Too goes 
to zero in the case with gravity, but in the case without 
gravity it goes to a negative constant. 

In the case without gravity, we can determine easily the 
value of the critical value of a, which turns out to be Oc = 
&^/(4c). But in the case of existence of gravity, just like 
the discussion in Ref. [14], the critical value of a is not Oc 
but a smaller effective critical value a, . The reason why 
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FIG. 6. The profiles of the background scalar field 4>- The parameters are set to & = 2, c = 1, and a = 0, 0.8, 0.836. 



< ttc can be interpreted as follows: When a = a» the 
thick brane will coexist in three phases, in other words, 
the solution of the background scalar (j) has the form of a 
double kink, which means (?i>'(0) = 0. Then from Eq. (23) 
with the conditions (24) we find this leads to the result 
that the bulk cosmological constant A = 0. So when 
y ^ oo, A'^ioo) = -Kim(a*)/3 (Knin(a) = V{(j){oo)) is 
the global minima of the scalar potential with fixed a. 
Because of the existence of gravity, i.e., A'^{oo) > 0, we 
have Knin(a*) < 0. Thus, noting that Vmin(a) increases 
with a and Vinin(ac) = 0, we reach the conclusion a* < ac- 
We also find just like the case without gravity when 
a = a* the solution of (j) is not unique. For 6 = 2 and 
c = 1, the critical value of a is a* = 0.837. 



by performing the coordinate transformation 

dz = e"^^y^dy. (27) 

The Dirac action of a massless 5D spin 1/2 fermion cou- 
pled to the scalar and gravity is 

where F^^ = {e~^^^,e~^j^) are the curved space 
gamma matrices, com are the spin connection [83-85] 
and their nonvanishing components are [■^)3, 73] 



B. Fermions on the brane with gravity 

The localization of fermions on the brane with gravity 
is very similar to the case without gravity. But in order 
to get the mass-independent potential for KK modes of 
fermions, we will follow Refs. [3, 5] to change the metric 
(20) to a conformally flat one 



'2A 



[rj^jjdx^dx^ + dz^ 



(26) 



where the prime denotes the derivative with respect to 
z. Using the conformal metric (26), the resulting Dirac 
equation is 



* 0. 



(29) 



7^a^ + 7'(9. + 2yl')-'/e^ 
With the usual KK and chiral decomposition 

■^[X,Z) = (i'l.n{x)hn{z)+1pKn[x)f^n[z)y{iG) 
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FIG. 7. The profiles of the warp factor e^^. The parameters are set to 6 = 2, c = 1, and a = 0, 0.8, 0.836. 





FIG. 8. The profiles of the matter energy density Too. The parameters are set to 6 = 2, c = 1, and a = 0, 0.8, 0.836. 



we can obtain the coupled equations for /l„ and /r„ [G.3] : 

[d, + 2A' + 7^e^F{(t>,x)]hn{z) = ™„/r„(^), (31) 
[d, + 2A' - 77e^F(0, x)] fnniz) = -m„/L„(z), (32) 

where fhn and /r„ satisfy the following orthonormality 
conditions: 

/oo 
e^^ hmfKndz = (5to„(5lr- (33) 
-oo 

Defining /l = /lc"^'^ and /r = /rb"^^, we obtain the 
Schrodinger-like equations for the new functions /l and 
/r [63]: 

[^dl + {z)\ /l„ = mi /l„ , (34a) 
{^dl + V^{z)\ fnn = mlfnn, (34b) 

where the effective potentials are given by 

Vi^iz) = 77^62^02 „ ^^A^, _ r,^e^A\ (35a) 
Vn{z) = ri^e^'^cf^ + r^e^4>' + r)4>e^A'. (35b) 

The equations of motion of the background scalar field 
and A in the z coordinate with metric (26) can be written 



as 

^" = _3^y + e2^^KM, (36a) 
d4> 

A" = A'2 - \<p'\ (36b) 

A'2 = V-ie2A(F(0)+A). (36c) 
6 3 

Before solving numerically the Schrodinger equa- 
tions (34), let us analyze the behavior of Vi^{z) and Vyi{z) 
at z ~ and z — > ±oo. 

At z 0, A(0) A'{Q) = </-(0) = 0, so Fl(0) = 
— ?70'(O) and Vr,(0) = rj4)'{jS). From the numeric solu- 
tions of (j){z) in Fig. 6, we can find (/)'(0) > 0. Then for 
77 > we reach the conclusion Vl(0) < and T^(0) > 0. 
When z — > ±00, (/)'(±oo) = 0. By the analysis to Eq. 
(36b), we see that A{±oo) oc — ln(±z). So ^'(±00) cx 
Tl/zU-^oo ^ and e2^(±°°) = l/^^U^oo ^ 0. There- 
fore, Vl(z) and Vr (z) fall off to zero at infinity. 

Equations (36) have been solved numerically with the 
boundary conditions (24) for A and 0. The profiles of the 
potentials (35) are shown in Fig. 9 for different values of 
a and fixed 6 = 2, c = 1, and 77 = 2. From the potential 
of the left-handed fermion shown in Fig. 9, we find that, 
with an increase of the value of a, both the height and 
width of the potential increase and there will be a double 
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FIG. 9. The profiles of tfie potentials for left- and right-handed fermions at the 2 coordinate. The three lines correspond to 
a — 0, 0.8, 0.836, respectively, and with parameters b = 2 and c = 1. 



well. The minima of the potential are always lower than 
zero, so the left-handed fermion may have a zero mode. 
The zero mode is plotted in Fig. 10 for different values of 
a. One can see that the width of the zero mode increases 
with the value of a. For the right-handed part, there 
is also a potential well. As the value of a is increasing, 
both the height and depth of the potential Vyi{z) are 
increasing. However, there is a difference from the left- 
handed part, i.e., the minimum of Vr is always positive. 
So the right-handed fermion will not have a bound zero 
mode. 

Figure 9 shows that the potentials of left- and right- 
handed fermions are indeed the modified volcano-type 
ones. Hence, there will only exist continuous spectrum 
for the KK modes of fermions. The volcano-type po- 
tential implies the existence of resonant or metastable 
states of fermions which can tunnel from the brane to 
the bulk [7f-73]. In order to get the solution of mas- 
sive KK modes /l,r(^), we need two initial conditions 
for Eqs. (34) at z = 0, which can be set as 



/(0) = ci,/'(0)=0 
for the even parity part and 

/(O) = 0, /'(O) = C2 



(37) 



(38) 



for the odd parity part, where Ci and C2 are arbitrary 
constants. Here we choose ci = 1 and C2 = 1. Since 
Eqs. (34) are Schrodinger-like we can interpret |/l,r(^)P, 
after normalizing fh,R{z), as the probability for finding 
the massive KK modes on the brane. But the massive 
modes cannot be normalized because they are oscillating 
when far away from the bane along the extra dimension. 
Ref. [72] proposed a function 



PL,R(m2) 



III \k,RM'dz 



(39) 



as the relative probability for finding the massive KK 
modes with mass square within a narrow range 
—Zf, < z < Z}, around the location of a brane correspond- 



ing to a larger interval 



< z < z„ 



where 2zb 



Zmax = 8zh here. So for KK modes with ^ W^'is'^j 
the values of Pl,r(™^) will be about 1/8 because /l.r 
can be approximately taken as plane waves. 

As examples, we have plotted the profiles of Pl,r(?ti^) 
in Figs. 11 and 12 corresponding to a = 0.8 and a = 
0.836, respectively, with b = 2 and c = 1. From the two 
figures we can see that the number of resonances is the 
same for both the left- and right-handed fermions and 
increases with the value of a. In these figures each peak 
corresponds to a resonant state. We can estimate the 
lifetime t of every resonance as r ~ where F = Sm 
is the width of half of the height of the corresponding 
peak [86]. For the case of a = there is no resonance 
found. 

The mass m, width F, and lifetime r of the resonances 
with different values of a arc listed in Table HI. From 
the Table we can see that the masses of left- and right- 
handed fermion resonances are almost the same, so the 
formation of the massive Dirac fermions can be realized 
[71, 72]. Table HI, Fig. 11, and Fig. 12 also show that 
the increase of the value of a will increase the number of 
resonances. 



IV. CONCLUSION 

In this article, we have investigated the effects of the 
variation of the mass parameter of a on the thick branes. 
We use a real scalar field, which has a potential of the cj)^ 
model, as the background field of the thick branes. 

In this model the critical value of a of the phase transi- 
tion on the brane with gravity will be lower than the one 
on the brane without gravity, and this result is consis- 
tent with the conclusion in Ref. [14]. The discrete bound 
modes of fermions in the case without gravity turn into 
quasilocalizcd resonant states in the case of with gravity 
[59]. 

We find when the critical value of a is reached the 
solution of the background scalar field is not unique, and 
it has a double-kink form. This happens in both cases 
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a 


Parity and Height of Vl,r 


n 




m 


r 


r 


0.836 


Left-handed 

y^max ^ 2.315306 


1 


0.307717 


0.554722 


0.00001623 


61614.2945 


2 


0.95514 


0.97731 


0.00213339 


468.7375 


3 


1.74224 


1.31994 


0.0260116 


38.4444 


Right-handed 
V^''''' = 2.396691 


1 


0.307716 


0.554722 


0.00001623 


61614.2945 


2 


0.95513 


0.97731 


0.00209757 


476.7421 


3 


1.73979 


1.31901 


0.0260411 


38.4008 


0.8 


Left-handed 

y^max ^ 2.23869 


1 


0.890670 


0.943753 


0.002798 


357.3981 


2 


2.02992 


1.42475 


0.24784 


4.0349 


Right-handed 
V^^'^'' = 2.324305 


1 


0.890670 


0.943753 


0.002850 


350.877 


2 


2.02980 


1.42471 


0.2928 


3.4153 



TABLE III. The mass, width, and lifetime of resonances of the left- and right-handed fermions. For the case of o = no 
resonance exists. The parameters are & = 2, c = 1, and r] — 2. V{^^ is the maximum of the height of the potential Vl,k, and 
n is the order of resonant states corresponding from small to large. 
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FIG. 10. The profiles of the zero mode for the left-handed 
fermion at the z coordinate. The three lines correspond to 
a = 0, 0.8, 0.836, respectively. Other parameters are set to 
h = 2 and c = 1. 



with and without gravity. That means at the critical 
value of a the width of the thick brane is not fixed, and 



with a disordered phase between two ordered phases [14]. 
Our other main result is that the number of the bound 
states (in the case without gravity) or the resonant states 
(in the case with gravity) increases with the value of a. 
This means that the branes with a big value of a would 
trap fermions more efficiently. 
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